Interaction-induced harmonic frequency mixing in quantum dots 
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We show that harmonic frequency mixing in quantum dots coupled to two leads under the influence 
of time-dependent voltages of different frequency is dominated by interaction effects. This offers a 
unique and direct spectroscopic tool to access correlations, and holds promise for efficient frequency 
mixing in nano-devices. Explicit results are provided for an Anderson dot and for a molecular level 
with phonon-mediated interactions. 
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The nonlinear mixing of two signals with different fre- 
quencies is a widespread and important concept used in 
many areas of physics. Frequency mixers have been ex- 
tensively used for a long time in microelectronics [3], 
where a diode provides the required nonlincarity, and 
in nonlinear optics (three- and four- wave mixing) 
Submillimcter-wave heterodyne signal detection at the 
quantum noise limit based on the nonlinear current- 
voltage characteristics of a superconducting tunnel junc- 
tion (ST J), is ubiquitous and has generated an enormous 
boost in radioastronomy. Apart from STJs, however, fre- 
quency mixing in quantum-coherent mesoscopic or nano- 
devices has been experimentally realized only for a three- 
terminal ballistic junction [J] and for the single-electron 
transistor [B], [|| , the latter possessing a broad and tunable 
range of frequencies and bandwidths. 

In this work, we give the quantum theory of frequency 
mixing in quantum dots - serving as prototype examples 
for nano-devices - in the presence of electron-electron 
(e-e) or electron-phonon (e-ph) interactions. Mesoscopic 
systems driven by ac fields are of major interest in con- 
densed matter physics @, i, @, EE EI El El EI El, El, 
E3, E3 > but a quantum theory of mixing has only been 
formulated for STJs @. We consider a multi-tone setup, 
where the dot is attached to source and drain contacts 
( "leads" ) with time-dependent (ac) voltages of different 
drive frequencies u>l /r (l9l | , and an additional gate is ca- 
pacitivcly coupled to the dot. A quantity of main interest 
is the time-dependent current I(t) through the device, 
expanded in a Fourier series as 



J(i) = Re J2 
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with frequencies uj nm = nui^ + "mujR. By determin- 
ing the complex- valued I nm , we derive general condi- 
tions under which I(t) exhibits harmonic frequency mix- 
ing. A mixing current can arise only if there exists at 
least one pair (n > 0, m ^ 0) with mixing amplitude 



In 



I-n -m 7^ (the- star denotes complex 



conjugation). In order to establish the importance of in- 
teraction physics on frequency mixing, we explicitly com- 
pute the J nm for two important models, namely (i) for 
on-site e-e interaction U in a spin-degenerate single-level 



dot (Anderson impurity model), and (ii) for a spinlcss 
level with e-ph coupling to an Einstein phonon, describ- 
ing transport through vibrating molecules [IE [2l[ . 

Before turning to derivations, let us briefly summarize 
our main findings in simple qualitative terms. Taking 
the standard wide-band limit (WBL) for the leads, we 
find that no frequency mixing is possible in the absence 
of interactions. This is a striking and unexpected result 
that we shall discuss in some detail. Once interactions 
are present, however, one will generally encounter mix- 
ing, J nm ^ (23 | . Therefore harmonic frequency mixing 
is dominated by interactions, and thus provides a highly 
sensitive novel spectroscopic tool of correlations. Impor- 
tantly, we show that even in the linear response regime 
(dc bias voltage V — ► 0) frequency mixing occurs in gen- 
eral, where the nonlincarity required for mixing is now 
generated by interactions. This implies that one may be 
able to operate such a quantum-coherent nano-device as 
frequency mixer with reduced dissipation and shot noise. 

We here take the left/right (a = L/R = ±) con- 
tact's single-particle energies as time-dependent, 6&, a — ► 
£fc. Q ± eV/2 + V£ c cos(w Q i), with dimensionlcss drive am- 
plitudes a a = V£°/u) a (we put h = ks = 1). The 
time-dependent current I{t) = [/l(^) — In(t)\/2 is then 
expressed in terms of A. a (t) = a a sin(uj a t) and Fermi 
functions / Q (w) = f(uj — aeV/2) describing the respec- 
tive contact, and the exact but unknown retarded (lesser) 
Green's functions (GF) G r (t,t') (G<(t,t')) of the dot gj. 
A time-dependent gate voltage then produces a similar 
term A g (£). For simplicity, we assume that the tunnel- 
ing amplitudes connecting the dot to the leads are not 
modulated by time-dependent voltages, and correspond 
to energy-independent hybridization matrices r Q in the 
space of spin and dot level indices. By a gauge-invariant 
generalization of Ref. Q, the time-dependent currents 
I a {t) for an interacting multi- level quantum dot can be 
written in the form 



I a (t) = -2eImTr[r Q G < (t,t)] 

- 2eImTr J dt' J ^/^^e'^'-''' 

x e -ifc(*)-*«(t')]r a G r (t,t'), 
where S a (t) = A a (t) - A g (t). 
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The formula ([2]) is manifestly gauge-invariant: adding 
an arbitrary time-dependent voltage to all (source, drain, 
and gate) electrodes leaves I a (t) unaffected. How- 
ever, Eq. ^ per se does not include displacement cur- 
rents [23f , which are generally necessary for conserva- 
tion of total charge in time-dependent quantum trans- 
port [HI, [IB, H3, Eli ■ Moreover, they generate a driving- 
induced energy renormalization of the dot states [ill , Il6| | , 
which has to be determined self-consistcntly. On gen- 
eral grounds, Fransson [lH ] has shown that for arbitrary 
drive frequencies, the sclf-consistently evaluated displace- 
ment currents can be separated into independent (left 
and right) parts. Thus, they cannot contribute to the 
mixing amplitudes JT]), and we shall disregard them from 
now on. Having established gauge invariance, we also put 
A s (t) = in what follows. 

Using Eq. @ , it can now be shown that in the absence 
of interactions, frequency mixing is suppressed by the 
parametrically small ratio Tr(T L / R )/D, where D is the 
electronic bandwidth in the leads. Therefore, in the wide 
band limit for the leads, no frequency mixing occurs for 
a nonintcr acting quantum dot, J nm = 0, regardless of 
how complicated its level structure may be. Technically, 
this follows because in the WBL the noninteracting GF 
G r {t,t') is independent of the ac drive, while Gq (t,t) 
depends on the left/right drive in an additive way only 
Q. One therefore encounters only Fourier components 
Inm ¥= for n = or m = in Eq. (J). The WBL for the 
leads is known to provide an excellent description when 
transport is dominated by states close to the Fermi level 
where T Q has only weak energy dependence. Note that 
the same conclusion holds true for [/l(£) + I R (t)]/2, as 
the "no mixing theorem" applies to II and I R separately. 
Below we focus on the J nm extracted from I(t) = [/^(t) — 
I R {t)]/2, and adopt the WBL, where J nm ^ can only 
be caused by interactions. Moreover, from now on we 
assume that only a single dot level is relevant and take 
symmetric hybridization Tl = F R = T/2. The respective 
generalizations are straightforward but do not yield new 
physics. The absence of mixing in the noninteracting 
limit may come as a surprise, since even in the WBL the 
IV curve of an undriven single-level dot is nonlinear. For 
instance, for a noninteracting resonant level at T = 0, one 
finds I(V) = ^tan~ 1 (el//2r). While this nonlinearity 
allows to generate higher harmonics and rectification |l4 , 
[lH, it does not create finite mixing amplitudes in Eq. (fTj). 
The latter arise from an effective "cross-talk" between the 
source and drain electrodes, which can only be mediated 
by interactions. 

The fact that interactions induce mixing can be 
demonstrated on the simplest possible level by first- 
order perturbation theory in U for the Anderson dot. 
The retarded self energy is then given by T, r (t,t') = 
Un(t)S(t — t'), with the time-dependent dot occupation 
n(t) = —iG^(t,t). Expanding in Bessel functions, a 
closed expression for n{t) in terms of the free GF and 
Fermi functions follows. Notably, in the driven case, the 
dot occupation and hence S r become time-dependent, 



while for the time-independent case, there is only a rigid 
shift of the dot level. The GF correction AG r = G r i: r G r 
now generates mixing, and after some algebra, we obtain 
the mixing coefficients to first order in U (n, m ^= 0), 

oo 

I nm = iU 2_j Jk+n(aL)Jk(aL)Ji+m(a>R)Ji(aii) (3) 

kl — — oo 

x [ F a(( l + m)ujR, loj R )Ff(-koj L , -kuj L - muj R ) 
- Fa((k + n)uj L ,kuj L )F^(-lu R ,-lu R -nuj L ) , 
with the auxiliary functions 

f$V )Wa )=r J ^.f R/L ( e )G r (e + u; 1 )G r /a (e + cj 2 ), 

which can be evaluated in closed (but lengthy) form. 
They obey the symmetry relations F r (uii,u>2) = 
F r (u>2,uii) and F a (u)i, uj 2 ) = F*(uj 2 ,uji). Elementary in- 
spection of Eq. ((3|) shows that for V — 0, I nm — > — I mn 
for u) R as expected when exchanging source/drain 

contacts. Eq. ([3]) nicely illustrates the basic mechanism: 
for U 0, terms containing lol and lo r appear in a multi- 
plicative way, showing that mixing is due to "cross-talk" 
of the leads. For zero dc bias (V = 0) and eo = 0, Eq. © 
and the above symmetry relations for F r / a imply J nm = 
for all even n + m. Frequency mixing in this particle- hole 
symmetric limit thus disappears, say, at the difference 
frequency cj^-i = lol — Wjj, but is still present at £J2,-i- 
However, once V ^ or eo ^ 0, one always finds mixing. 
To see this, it is instructive to consider very small drive 
amplitudes V£ C R in Eq. ([3]). Upon expanding the Bessel 
functions, the mixing amplitudes J nm cx a^a R are non- 
zero only for n,m = ±1, and the double sum ([3]) receives 
contributions only from k = 0, — n and / = 0, — m. Since 
n + m is now always even, there is no mixing unless V =/= 
or eo ^ 0. This simple calculation already demonstrates 
that mixing is generated by interactions. 

To study interaction effects beyond lowest order, we 
next consider the sequential tunneling regime with T > 
r, where a master equation approach applies. Here one 
evaluates the dynamics of the occupation probabilities 
P s (t) [with < P s (t) < 1 and J2t=i p s(t) = 1] for the 
four possible dot configurations s with energy e s : s = 1 
denotes the empty dot (ei = 0), s = 2,3 the singly- 
occupied dot with spin up/down (e 2 .3 = eo), and s — 
4 the doubly occupied one (€4 = 2eo + U). Following 
standard steps [lOj . the master equation reads P(t) = 
T. s -, a =L/R\ K sAi)PAt) - K« s (t)P s (t)] with transition 
rates (s ^ s' , where K° = - J2s'^s K Z>) 

OO 

K«,{t) = T a Rc J2 e^ t J,(a Q )J fe+? (a Q ) (4) 

k,q— — 00 

x \^N ±g± {±[e s -e s ,] + kuj a - aeV/2), 
± 

where N + = J2a=i 1 K s Mo-l s ')| 2 i N- follows by replacing 
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FIG. 1: Mixing amplitudes Ji,±i in the sequential tunneling 
regime, in units of 2eT/h for V = e = 0,lj l = 20F,V£ C = 
Vr c = 50r and T = 5r, as a function of (a) lu r /lu l for 
C7 = 100r, and (b) of U/T for lj r = 50r. 



the dot fermion operator d\ —> d a , and 

„,n i , -D ( 1 ie 

fl± e = /(e) ± - In — T -Re V r + 5-= 

7T Z7TJ 7T V 2 Z7TJ 



, (5) 



where -0 is the digamma function, implying logarithmic 
divergencies }24| . After Fourier expansion as in Eq. |Tj), 
the master equation leads to an algebraic equation. Nu- 
merical solution obtains the Fourier coefficients P s .nm 
and hence the J nm from I a (t) = J2 S ^ S ' ^as'K^,(t)P a >(t), 
where ss i = ±1 for s ^ s'. These equations allow 
to reproduce known results for the sequential tunneling 
current under a dc bias and for the ac driven case with 
&r = ujl QUI- F° r U = 0, in accordance with our dis- 
cussion above, no frequency mixing is found from the 
master equation, while J nm ^ for U ^ 0. To ensure 
consistency, we have checked that for small U, the master 
equation reproduces the perturbative result ([3]) taken at 
high T or large V. 

Figure [T] shows the mixing amplitudes Ji,±i for V = 
eq = as a function of U/T and lo^I^l- Note that for 
uj r = lul, we correctly find J\\ = as enforced by the 
V = symmetry I nm = —I mn under exchange of ui R and 
u>L. At this point, ReJi._i also vanishes, while ImJ! j 7^ 
does not generate current, 7(f) oc sin[(cjR — w^jt] = 0. 
Remarkably, the mixing amplitudes display characteristic 
features (peaks or steps) at certain ratios Ufl/wi, cp. for 
ur/ujl = 1/3 in Fig. [Tfa). Such features can be rational- 
ized in simple terms by noting that the ac voltages corre- 
spond to photon- assisted side peaks in the dot's density 
of states, located at energies eo ±w nm and 2cq + U ±oj nm 
with arbitrary integers n, m. Once one of those energies 
hits the Fermi level (which is located at eo = in Fig.[Tj), 
transition amplitudes are resonantly enhanced, and mix- 
ing becomes particularly efficient. To give an example, 
such resonances occur for ljr/lul = (U /uji, — n)/m, and in 
Fig. [IJa), where U/ojl = 5, the feature at u) R juj L = 1/3 



corresponds to (n, m) — (4, 3). Which of these commen- 
surability features (indexed by n, m) will actually show 
up in the mixing current is primarily determined by the 
drive amplitudes V£J R . For the rather large V£J R taken 
in Fig. [TJ large (n, m) have to be taken into account, 
resulting in the rather complicated dependence of the 
mixing amplitudes on u r /ujl observed for lu r /lul < 1/3. 
The data in Fig. [lja) show that mixing disappears in 
the limit col —> 00. We also observe that mixing dis- 
appears in the opposite limit u> R /u>i — > 00, where one 
can effectively average over the fast lo r oscillations and 
ends up with a monochromatic situation again. Finally, 
additional calculations (not shown) reveal that like in 
the small-?/ case, there is no mixing (J nm = 0) at the 
particle-hole symmetric point, eo = — U /2 and V = 0, for 
all even n + m. Before turning to the case of e-ph interac- 
tions, we briefly comment on the Kondo regime, realized 
at low temperatures for U T and eo ~ —U/2, where 
one can map the Anderson dot to a spin-1/2 impurity 
problem. At the special Toulouse point, an exact solu- 
tion of the mixing problem can then be obtained from 
Ref. [l3||. This solution will be discussed elsewhere, but 
by taking T = V = 0, mixing is seen to disappear, in ac- 
cordance with the Fermi liquid nature of the Kondo fixed 
point. At finite T, V, however, we again find nonzero 
mixing amplitudes. 

Next we discuss a spinless level (fermion operator d), 
with e-ph coupling A to a phonon mode Q = b + b* of 
frequency VL (boson operator b). The dot is described by 
-ffdot = [eo + XQ]d^d + ilb' !f b, and connected to leads as for 
the Anderson dot. We have computed the mixing ampli- 
tudes Jnm = Inm + I- n ,-m by second-order perturbation 
theory in A, with the result I nm = [Inm ~ ^nl/^i where 

C = »(Ar/27r) 2 Ji +n (a a ) Ji(a a ) Jfc +TO (a_ a ) J k {a~ a ) 
k,l 

x J du> J def a (w)f- a {e)Dl(u) - e-lwa-ku-a) 

x G r (e + (k + ™)w_ a )Go(e + kuj^ a )G r {cj - lto a ) 
x Gq(oj - lui a + mu)- a ), (6) 

where D r (u) = 2lj/[{lj + i0+) 2 - O 2 ]. From these ex- 
pressions, one shows that for V = eo = 0, again J nm 
only for even n + m. For f2 ^> T, the high phonon fre- 
quency allows to approximate the boson propagator by 
a constant, and one essentially comes back to the Ander- 
son dot result ([3]) where mixing was established above. 
In Fig. [51 we show the mixing amplitudes for Q = T and 
small drive amplitudes in Eq. ([6]). 

Finally, we comment on how frequency mixing could 
be verified experimentally, and on consequences for ap- 
plications. Working quantum-dot frequency mixers 0, [|| 
operate in the adiabatic regime at rather low lu R /l (no 
photon-assisted effects) , probably outside the applicabil- 
ity range of our theory [19(. The measured dot capaci- 
tance C = 550 aF 0] puts U = e 2 /2C around 35 GHz, 
while r = 2-k/RC w 13 GHz follows from the measured 
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FIG. 2: Mixing amplitudes J nm for phonon-mediated in- 
teractions and small drive amplitudes V£J R , see Eq. (JS| ex- 
panded to order V£ c Vr°. Parameters: T = e = 0, ft = 
r, V — 0.8F, u>l = 2r. The J nm are given in units of 
\ 2 V£ c V£ c /(2ivry i . 

resistance R = 850 kfi, implying U/T 2.7. By com- 
paring to Fig. [1] a large mixing signal is expected for 
U/T w 50, which may be achieved for the same C by 
increasing i? to R rts 23 MX2, leading to T = 500 MHz. 
Choosing lo r / l in the GHz range, see Fig. [TJ should then 



put the device into an experimentally accessible regime, 
where our theory applies and predicts strong frequency 
mixing. In fact, one may be able to generate THz waves 
by such a mixer. In the present work, we have shown 
that mixing is possible for V = 0, where the dc current 
vanishes and hence a very small noise level is expected. 
Existing experiments [j| involved finite dc current, caus- 
ing large noise levels and dissipation effects in the de- 
tection electronics. Both these problems can be consid- 
erably reduced by using a correlated dot near V = 0, 
where interactions provide the required nonlinearity. 

To conclude, the theory of harmonic frequency mix- 
ing in interacting quantum dots has been given. In the 
wide-band limit for the leads, a "no mixing theorem" can 
be established, stating that mixing requires the presence 
of interactions. For both e-e and e-ph interactions, we 
have then shown that mixing is indeed generated, and 
provided detailed quantitative predictions for the mix- 
ing amplitudes. We hope that these findings stimulate 
experiments and further theoretical developments. 
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